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Abstract 

We demonstrate that the one-dimensional (ID) transport channels that appear in the gap 
when graphene nanoroads are embedded in boron-nitride (BN) sheets are more robust when 
they are inserted at AB/BA grain boundaries. Our conclusions are based on ab-initio electronic 
structure calculations for a variety of different crystal orientations and bonding arrangements at 
the BN/C interfaces. This property is related to the valley-Hall conductivity present in the BN 
band structure and to the topologically protected kink states that appear in continuum Dirac 
models with position dependent masses. 

Metallic transport channels appear at the edges, surfaces, and interfaces of two and three di- 
mensional bulk insulators when a bulk topological index changes value as the interface region 
is crossed.- - — This property can provide transport channels in otherwise insulating materials. The 
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metallic states possess an internal structure related to their sense of propagation which leads to 
special transport properties including zero bend resistance at sharp turns in the current propagation 
trajectory, pseudospin memory, and suppressed backscattering. 7 



Same topology 

Zigzag 



o b Opposite topology 

Zigzag 



Armchair 



Armchair 



H - (BN) n - C m - (BN) n - H H - (BN) n - C m - (NB) n - H 

Figure 1: (Color online) Schematic plot of two types of zigzag and armchair BNC nanoroads. 
The integer n is the number of atoms in each BN partition whereas m is the number of C atoms 
across the nanoroad. The outermost BN atoms are passivated with hydrogen. Left Panel: The same 
topology arrangement in which the sign of the valley Hall effect is the same in both BN sheets. In 
the absence of C atoms the BN fragments can be joined seamlessly. Right Panel: The opposite 
topology arrangement in which the sign of the valley Hall effect is opposite in the two BN sheets. 
In the absence of C atoms it is necessary to form a nearest-neighbor bond between atoms of the 
same species- either N— N and/or B— B. 



In this Letter we show that two-dimensional hybrid structures consisting of graphene nanoroads— 
(See [[]) embedded in hexagonal boron nitride (BN) sheets can be an attractive host for topologi- 
cally assisted one-dimensional (ID) transport channels. Our study is motivated by the observa- 
tion that the 7T-bands of BN are similar to those of a graphene sheet, except that the honeycomb 
sublattices have different ^-electron site energies. (The electrostatic potential is more attractive 
on the higher-Z N atom sites). Because valence band states near the K and K' Brillouin-zone 
corners are strongly localized on N sites whereas those elsewhere in the Brillouin-zone are di- 
vided more evenly, the bands of BN contribute with valley Hall conductivities and associated 
Berry curvatures of opposite signs near K and K' whose continuum model expression is given 
by Q(q) = TV (3a 2 At 2 ) / ^2 (A 2 + 3q 2 a 2 t 2 ) 3 ^ 2 ^ where q is the momentum measured from a Dirac 
point, x z = ±1 is the valley, a is the lattice constant of the honeycomb lattice and A is the band 
gap. 14 If the B-site vs. N-site potential difference was small, these Berry curvatures would be 
strongly concentrated near the BZ corners [see [21 . Indeed, in the graphene case it is known that 
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Figure 2: Berry curvatures obtained from a tight-binding model of the honeycomb lattice with 
different band gap sizes A = 2U where U is the staggering site potential. The Berry curvatures 
are more sharply peaked near the Brillouin zone corners K and K' when the band gap is small and 
they spread out more when the gap is larger. The rightmost panel shows the deviation of the valley 
resolved Chern number from the ideal 1/2 value of the continuum model obtained integrating the 
Berry curvature in one equilateral triangle corresponding to one half of the primitive cell associated 
to valley K given by Ck = (2tt) _1 j K d 2 k£l(\<L). The dotted vertical line at 2U jt = 2 indicates the 
valley resolved Chern number obtained for the tight-binding hamiltonian we used to approximate 
the BN bands. 



electronic properties in systems with sublattice- staggered potentials that are weak compared to 
the 7T-band width can be described using a two-dimensional massive Dirac equation. Under these 
circumstances Dirac-equation continuum models are valid and predict one-dimensional (ID) local- 
ized states along lines where the sublattice-staggered potential (the Dirac-equation mass) changes 
sign. 3 ' 4 In the closely related bilayer 5 " 7 and multilayer 6 graphene cases, the sublattice staggered 
potentials are readily generated experimentally by applying an electric field across the layers and 
varied spatially by appropriately arranging external gate voltages. In the single-layer case, how- 
ever, it has not been obvious how the sublattice-staggered potentials could be realized, although 
gaps may be present with a lattice matched BN substrate 11 or when gate potential profiles are 
carefully correlated with strains. 12 This work is motivated by the idea that the difference between 
B-atom and N-atom potentials in BN provides the desired staggered potential. Of course the stag- 
gered potential is not weak in the case of a BN crystal, so that expectations based on massive Dirac 
continuum models must be checked by ab initio electronic structure calculations. 

The topological index that is inherent in BN bands can be viewed as a valley Hall effect since 
each valley separately supports half integer quantum Hall effects of opposite sign, or at least does in 
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the limit that the potential difference is small compared to the 7T-band width. The properties that we 
discuss below are therefore closely related to the simple momentum space Berry curvature pattern 
in the bands of BN sheets. Hybrid BN/graphene systems of the type we study occur naturally 
in patched sheets containing a mixture of atomically thin graphene and BN, 15 and have recently 
become a subject of great interest,-^*^ - — 2Sr21 as the basis of a possible strategy for controlling 
graphene system band gaps. In this work we show that the properties of the transport channels 
formed by graphene nanoroads in BN strongly depend on the valley Hall effect of the surrounding 
material. 

The simplest graphene nanoroad geometries have C/BN interfaces with crystallographic ori- 
entations along the zigzag and armchair directions. The unit cells considered in our ab initio 
calculation are illustrated in[IJ In each case a graphene nanoribbon is flanked on the left and right 
by BN sheets with B and N atoms either on the same or opposite sublattices. When the two dif- 
ferent BN fragments can be joined seamlessly preserving the crystalline B and N atom sequences 
upon carbon atom removal, as in the cases illustrated in the left panels of [Q the BN sheets on 
left and right have valley Hall conductivities of the same sign and we will say that they have the 
same topology. The right and center panels show junctions in which the BN sheets on opposite 
sides of the graphene nanoroad have opposite valley Hall conductivities. We will refer to these 
configurations as having opposite topologies. When the carbon atoms are removed, joining the 
BN fragments would in this case require B— B or N— N bond. As shown in[3l a difference in the 
topology of the BN sheet arrangement on the nanoroad shoulders invariably leads to a qualitative 
change in nanoroad quasi- ID band structure. 

The indices n and m in [3] represent respectively the number of BN and carbon atoms in the 
ID unit cell. The width of the BN fragment, specified by the index n, plays a minor role in 
mid-gap state properties. We begin by discussing results for zero, one, and three hexagon width 
zigzag nanoroads presented in the upper panels of [3] The zero width cases illustrates the electronic 
structure of seamless BN nanoribbons with In = 80 atoms of BN, and therefore a width of 86 
nm. The large band gap of 4 eV is comparable to that of bulk BN. There is a notable lack of 
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particle-hole symmetry in the band structure, which is partly due to the difference between B 
and N site potentials mentioned previously but cannot be captured by a nearest neighbor tight- 
binding model. We are mainly interested in the mid-gap states which emerge when a graphene 
nanoroad is embedded in the BN sheet. As expected on the basis of the properties of pure graphene 
nanoribbons, we find that zigzag graphene nanoroads do have strongly reduced gaps because of 
the presence of zigzag edge states.— - — There is still however a topology sensitive feature; in the 
opposite topology case a mid-gap state appears that is peaked near the center of the nanoroad, 
rather than at its edges, with some penetration into the BN region. This state is indicated by blue 
shading in Figs. 3(d) to 3(i). We view this state as the lattice remnant of the mass-reversal (kink) 
state which would appear in the opposite topology case if the staggered potential was weak and the 
continuum model applied. Even in the absence of carbon atoms its appearance reduces the energy 
gap by approximately 2 eV. If the continuum model was accurate there would be opposite velocity 
kink states associated with the K and K' valleys and the gap would disappear entirely. 

Although gaps are absent for wide nanoroads in both like cases, the additional kink state present 
for opposite topologies closes the gap already at one hexagon width in the opposite topology case. 
For same topology configurations, the gaps are closed because the k x ~ states that are most 
strongly localized at the edge are shifted in opposite directions by coupling to the BN shoulders.— 
For inverted topology, since the same type of atoms N or B are attached at both edges of the 
graphene nanoroad, we see a shift of onsite energies at the edge sites in the same direction, upwards 
for N atom bonds and downwards for B atom bonds. In this case the states closing the band gap 
have kink- state character and are spread widely over the C atoms of the nanoroad, and with some 
penetration into the BN region. 7 

We now turn to armchair nanoroads. In ribbons this orientation does not support metallic 
edge states, 40 and instead yields gaps which scale as the inverse ribbon width. Continuum model 
considerations suggest that reversed topology would lead to vanishing gaps for all widths. 6 ' 7 The 
band structures in the lower panel of[3]demonstrate that there are gaps in the reversed topology case, 
but that the gaps are much smaller than in the same topology case. The magnitude of the small 



5 



avoided crossing gap between opposite velocity kink states at k x ~ shows the same oscillatory 
decline as a function of ribbon width that has been heavily studied in ribbons. -22>22 The gap reaches 
a value as large as 0.75 eV for armchair single-hexagon nanoroads. It is noteworthy that this 
avoided crossing is always smaller than the finite-size gap in either nanoribbons of the same width 
or same topology nanoroads, a fact that allows the kink channels to become the main current 
conducting path when the Femi level of the system is adequately shifted. The same topology and 
opposite topology gaps, A and 8, are plotted as a function of nanoroad width in the bottom right 
panel of [H 

Zigzag graphene nanoroads have bearded interfaces terminations when they have an odd num- 
ber of carbon atoms or non modulo 4 even integer number of C atoms in the unit cell.— £L These 
interfaces are not expected to be common in experimental systems because of stability consid- 
erations. A combination of zigzag and armchair fragments 35 are more abundant than bearded 
configurations with repeated dangling atoms. Bearded edges lead to more complicated band struc- 
tures that we discuss further in the supplementary information, but generally preserve the tendency 
toward greater conduction in the opposite topology case. 

To address the influence of topology on nanoroad transport, we consider a two-terminal geom- 
etry with a central scattering region containing nanoroads that are either straight, or curved with 
both armchair and zigzag segments [See0. In our model, the source and drain electrodes are semi- 
infinite zigzag graphene nanoroads. The central scattering region contains 96x60 lattice sites and 
the nanoroad width is fixed at 32 carbon atoms. The BNC nanoroad system is modeled using a 
nearest neighbor tight-binding Hamiltonian of the form 

H = - f E C 1 C J + E U Ach - £ Uacja, (1) 

(ij) i<EA ieB 

where cj (c ; ) is a ^-orbital creation (annihilation) operator for an electron at the site i, Ua and Ug 
are the BN ^-orbital site energies on the A and B subllices, and t = 2.6 eV is the nearest neighbor 
hopping energy for both graphene and h-BN. We introduce a staggered AB sublattice potential 



6 



to describe the site energy differences between B and N sublattices of h-BN, i.e. Ub = —U^ = 
XUq. Here, A='±' defines the crystal topology of the h-BN regions in the graphene nanoroad. Uq 
measures the strength of the staggered potential, and for the purposes of this illustrative calculation 
we have assumed it to be Uq = t. 

The temperature T = conductance from the left lead to the right lead can be calculated using 
the multi-probe Landauer-Biittiker formula: 41 

e 2 

G RL = -Tr[T R GT L G a ], (2) 
n 

where G r,a are the retarded and advanced Green functions of the central scattering region. The 
quantity T R / L is a linewidth function describing the coupling between left/right-lead and the central 
region and is obtained by calculating the self-energy Z r R > L of left/right-lead using a variant transfer 
matrix method. 4 -"^ 4 - 

|5] plots conductance Grl as a function of Fermi energy £ for the straight (zigzag-edged) and 
curved BNC nanoroads. In panel (a), the crystal topologies are the same in upper and lower h- 
BN regions. We observe that the conductance of the straight BNC nanoroad is quantized to be 
G=2n e 2 /h (n = 0, 1,2...). The conductance plotted here is per spin, so that the even conductance 
values can be attributed to the presence of ballistic states that are localized at both edges. The 
nanoribbon finite size effect for this nanoroad width produces a small band gap with vanishing 
conductance near the charge neutrality point. When the graphene nanoroad becomes curved in 
the central regime, disorder produces backscattering which supresses the conductance below the 
quantized value. 

The case of opposite BN shoulder topologies is illustrated in[5£b). In this case the conductance 
rule for straight BNC nanoroads is changed to G = {In — 1) e 1 /h (n = 1,2...) as shown in Ob). 
This change reflects the addition of a topologically supported kink state. For the curved reversed 
topology nanoribbon, the conductance does not decrease substantially in the region near the charge 
neutrality point where the only remaining channel is the kink state. Near the neutrality point the 
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conductances show a surprising plateau feature with minor fluctuations that may be suppressed 
further by increasing the width of the graphene nanoroad. This finding strongly reflects the robust- 
ness of the crystal topology induced ID topological conducting state against BNC nanoroad bends 
and is consistent with the finding of the zero bend resistance for ID topological confinement states 
in gated bilayer graphene. 7 

To visualize the ID topological conducting channel, we can evaluate the position dependence 
of its local density of states (LDOS) of the current originating from P th lead using 

Pp(r,e) = ^[G'T p G fl ] rr . (3) 

Here £ is the Fermi energy. In [61 we plot the LDOS distribution of the topologically assisted 
conducting channel from the left lead into a snake-shaped BNC nanoroad. Our calculation shows 
that it is localized across the graphene atomic sites with only a small portion on the interface 
between h-BN and graphene. The transport calculations for the straight and curved nanoroad 
geometries confirms that the kink states picture remains essentially valid even when both zigzag 
and armchair edges are simultaneously present. The almost ballistic transport is preserved in spite 
of irregular bonds at the graphene/BN interfaces. 

In summary, we have shown that graphene nanoroads embedded in BN sheets host kink-like 
transport channels when surrounded by opposite topology BN fragments. Transmission by the 
kink-states is topologically assisted and nearly perfect even when the graphene nanoroad is curved 
and has irregular edges consisting of combinations of both for zigzag and armchair edges. 

Calculation Details 

The self-consistent ab-initio calculations were carried out using the PAW-LDA parametrization 
of the exchange correlation functional in a plane wave basis code, as implemented in Quantum 
Espresso. 45 We used the 1.443 Angstrom for the separation between neighboring B, N and C 
atoms and 1.09 for B— H and N— H separations. These values are consistent with optimized lattice 
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constants obtained from total energy minimization. Relaxation to obtain a final crystal geometry 
with atomic forces below 0.01 eV / Angstrom had a negligible effect in modifying the interatomic 
distances and band structures. For our nanoroad calculations in ribbon geometries we used 30, 
1 and 1 k-points in the k x , k y and k z directions, and a plane wave energy cutoff of 60 Ry. The 
vacuum separation distances between supercell repeated nano ribbons was set to 20 in both y and 
z directions. 
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Figure 3: (Color online) Upper Panel: Quasi ID band structures of zigzag (ZZ) graphene 
nanoroads of different widths embedded in hexagonal boron nitride. The three band structures are 
for widths corresponding to zero, one and three complete hexagons or equivalently m = 0, 4 and 
12 carbon atoms in a unit cell. The BN width number n = 40 for all cases. We consider the same 
topology configuration, (BN)„— C m — (BN)„, the reversed topology configuration A, (BN)„— C m — 
(NB)„, in which nitrogens bind with carbon and reversed topology configuration B, (NB)„— C m — 
(BN)„, in which borons bind with carbon. States that are peaked in the middle of the nanoroads 
are represented by thick blue lines. Lower Panel: Band structures of armchair (AC) graphene 
nanoroads of different widths embedded in hexagonal boron nitride. The two band structures are 
for m = 0, 4 and 12 carbon atoms in a unit cell, corresponding to nanoroads that are 0, 1 and 3 
hexagons wide. Same and reversed topology configurations are considered. The states that are 
localized in the middle of the nanoroads are again indicated by thick blue lines. As in the case 
of graphene nanoribbons the gaps have an oscillatory behavior with period 12 (corresponding to 3 
hexagons), embedded in a smooth envelope that follows a mT l law. 
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Figure 4: Schematic plot of the two-terminal BNC nanoroad transport measurement that we model, 
(a) Straight BNC nanoroads; (b) Curved BNC nanoroads: the two leads are exactly the same as 
those in panel (a). The nanoroad will support kink-like states when the two BN regions have 
opposite topology. 




Figure 5: (Color online) Conductance calculations for straight zigzag-edged and curved graphene 
nanoroads embedded in the h-BN as a function of the Fermi energy e. Panel (a): conductance 
for graphene nanoroads embedded within h-BN sheets with the same topology. An insulating gap 
clearly arises near the charge neutrality point. Panel (b): graphene nanoroads embedded within 
h-BN sheets with different topologies. The conductance plateau near the charge neutrality region 
is due to a topologically supported edge state. Irregular and curved edges introduce a extremely 
small backscattering probability. 



14 



-6 





x (nm) 



6 



Figure 6: (Color online) Local density-of-states of a kink state in a curved carbon nanoroad em- 
bedded between opposite topology BN sheets, suggested here by the '+' and '— ' signs, calculated 
through a nearest neighbor tight-binding model as discussed in the main text. The edges consist of 
a mixture of zigzag and armchair fragments but the disorder does not degrade transmission. 
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